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ABSTRACT 

It has been proposed recently that the small size and long lifetime of 
ultracompact (UC) HII regions could be due to pressure confinement if the 
thermal pressure of the ambient gas is higher than previous estimates. We 
point out that confinement by thermal pressure alone implies emission measures 
in excess of observed values. We show that turbulent pressure, inferred from 
observed non-thermal velocities, is sufficient to confine UC HII regions and 
explain their longevity. We predict an anti-correlation between the size of UC 
HII regions and the velocity dispersion of the ambient neutral gas, and show 
that it is consistent with existing observations. 

Subject headings: ISM: HII regions — ISM: kinematics and dynamics — stars: 
format ion — t urbulence 



1. INTRODUCTION 

Ultracompact HII regions (UCHII's) are very small (< 0.2 pc) dense regions of 
ionized gas in molecular clouds first noted by Ryle & Downes (1967) and Dreher & Welch 
(1981). Deeply embedded, the ionized gas is seen only through free-free emission at radio 
wavelengths. The importance of this type of object for our understanding of the birth of 
massive stars is well recognized (Churchwell 1993; Welch 1993; Vogel 1994). A key unsettled 
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question is the discrepancy between the characteristic age (~ 10^ yr) of UC HIIs imphed 
by their ubiquity (Wood & Churchwell 1989, hereafter WC89; Kurtz, Churchwell & Wood 
1994, hereafter KCW94) and the very short dynamical hfetime (~ 10^ yr) estimated from 
the sound crossing time. This disparity imphes the presence of physical mechanisms which 
confine the UC Hlls. Proposed possibilities include: (1) ram pressure due to gravitational 
infall of ambient gas onto the HII region (cf. Reid et al 1981; Dreher & Welch 1981; Ho 
& Haschick 1986), (2) stellar wind bow shock/ram pressure confinement (cf. WC89; van 
Buren et al 1990), and (3) photo-evaporating disks (Vogel, Genzel & Palmer 1987; Welch 
1993; HoUenbach et al 1994) or photo-evaporating clumps (Lizano & Canto 1995; Williams 
et al 1996), or champagne- fiow/blisters (Tenorio-Tagle 1979; Yorke et al 1983; Forster et al 
1990). 



2. CONFINEMENT BY THERMAL PRESSURE 



De Pree et al (1995) (see also Garcia-Segura & Franco 1996 and Akeson & Carlstrom 
1996) have recently suggested that the compactness of the UC Hlls can be explained if 
the density and temperature in the ambient neutral gas are much higher than commonly 
assumed. Their basic argument is two-fold. First, the initial Stromgren radius Rs scales 
inversely with gas density (e.g. Spitzer 1978) as 

^'~HixP^^ -1.99X10 (^Q49^-J 405cm-3^ 

where 5* is the flux of ionizing photons, no = 2nH2 is the initial electron density in the 
ionized gas and [32 = 2.6 x 10^^"^ cm^ s^^ is the recombination coefficient. Thus, for a 
sufficiently high ambient molecular density, the initial Stromgren sphere can be very small. 
Second, a higher ambient gas temperature reduces the amount that the ionized gas expands 
before reaching pressure equilibrium with the ambient material. 

The emission measure of the ionized gas predicted by this simple scenario, however, 
is significantly larger than the values generally observed. This point can be verified easily 
from the calculations and discussions by De Pree et al (1995). For a sample of small UC 
HIIs, De Pree et al showed that an ambient molecular density of 10^ cm~^ together with a 
gas temperature of ~ 200 K would suffice to explain the compactness of the UC HIIs. They 
estimated an initial Stromgren radius Rs = 1.0 x 10~^ pc at such high density, corresponding 
to an initial emission measure EMq = 2n'^Rs = 8 x 10^^ cm~^ pc. Overpressured, the initial 
ionized sphere expands outward until it reaches pressure equilibrium with ambient neutral 
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gas. At a later time, when the ionized sphere has a radius R and electron density rig, the 
emission measure will drop to 

EM = 2nlR = {^fEMo. (2) 
R 

For R = 0.01 pc, the emission measure is thus expected to be ~ 8 x 10^ cm~^ pc. The 
observed average emission measures for spherical UC HII regions with sizes ~ 0.01 pc are 
mostly around a few xlO^ to a few xlO^ cm~^ pc in the WC89 and KCW94 surveys, 
lower than the predicted value by over an order of magnitude, although there clearly are 
some UC HIIs which have very high emission measures (cf. Turner & Matthews 1984; 
Garay et al 1993a). It seems unlikely that the assumption of optical thinness for the radio 
emission could have underestimated the emission measure by such a large factor. Therefore, 
although the gas density around some UC HIIs is indeed higher than 10^ cm'^ (the density 
originally adopted by WC89) as recent observations have revealed (cf. Huttemeister et al 
1993; Akeson & Carlstrom 1996; Plume et al 1996), a gas density of nj^j = 10^ cm~^ as 
adopted by De Pree et al (1995) seem too high for the material around typical spherical 
UC HIIs with R ~ 0.01 pc in WC89 and KCW94 surveys. Other physical processes such as 
swept-up shells or dust absorption of Lyman photons may modify the predicted emission 
measure somewhat, but it is not difficult to see that the former would further increase the 
emission measure, while the latter can only reduce the emission measure by a factor of 
(1 — xY^^, where x is the fraction of Lyman photons absorbed by dust grains (Franco et al 
1990; Churchwell 1993). For x = 0.9, the emission measure is reduced by a factor of 0.46 
(WC89). While not necessarily irrelevant, dust absorption alone is not likely to be able to 
account for the observed low emission measure for a significant fraction of UC HIIs. We 
conclude that confinement by thermal pressure requires UC HIIs to have higher emission 
measures than observed in most cases. 



3. TURBULENT PRESSURE CONFINEMENT 



3.1. Plausibility 



That the gas pressure in the interstellar medium is not limited to the thermal pressure 
alone is evident from emission line profiles and has recently been brought up in connection 
with UC HIIs by Garcia-Segura & Franco (1996). Ever since the discovery of molecular 
clouds, it has been known that the Doppler-broadened line profiles of all molecular tracers 
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indicate the presence of significant non-thermal motions, referred to as turbulence, down 
to the smallest scales probed even for low mass molecular structures (cf. Falgarone, Puget 
& Perault 1992). Since massive star forming regions are well-known to have stronger 
turbulence (cf. Plume et al 1996), it is likely that turbulence exists ubiquitously on the 
scales of UCHIIs. Regardless of whether the origin is hydrodynamic or hydromagnetic, the 
presence of turbulence implies a significant pressure in addition to the thermal pressure. 
Massive star forming regions, where UC HIIs are found, are particularly turbulent, often 
with velocity dispersions well in excess of a couple of kilometers per second (Cesaroni et al 
1991; hereafter CWKC; Plume et al 1996). The turbulent pressure implied is considerably 
larger than the thermal pressure. For example, the turbulent pressure P = nH2f^H2^v 
corresponding to a velocity dispersion of = 2 km is equivalent to a thermal pressure 
at an effective kinetic temperature of Tf,ff ~ 10^ K. In fact, if massive star forming cores 
where UC HIIs reside are indeed close to virial equilibrium (CWKC), the implied molecular 
gas pressure ought to be comparable to the ram pressure due to free-fall collapse. Therefore, 
as long as turbulence exists down to the scales of UCHIIs, it seems safe to conclude that 
turbulent pressure could provide significant opposition to the expansion of the ionized gas, 
just as would the ram pressure due to gravitational collapse of the molecular envelope (Reid 
et al 1981). 

The validity of the above hypothesis can be clarified by further considering the two 
conceptual phases of UC HII development (Dyson & Williams 1980; WC89). The first 
phase, i.e., the establishment of an initial Stromgren sphere, is not affected much by the 
turbulence, so the general results and conclusions in Dyson & Williams (1980) and WC89 
remain valid, and the Stromgren sphere with radius determined by Equation (1) will form in 
a matter of a few years [|. The temperature of the newly ionized gas jumps to Th+ ~ 10^ K 
and the number density of the particles in the gas increases by roughly a factor of 4 due 
to dissociation and ionization. This, together with the ram pressure due to stellar winds, 
causes the ionized gas to expand outward at its sound velocity until reaching approximate 
pressure equilibrium with the ambient molecular gas. Instead of completely neglecting the 
pressure due to stellar winds and turbulence in the ionized gas (Dyson & Williams 1980; 
De Pree et al 1995), we introduce a factor ^ (> 1) to partially take them into account (It 
can be argued that the ram pressure due to stellar winds may be reflected by an enhanced 
actual electron density in a swept-up shell). The final average electron density rig in the 
HII region can then be estimated from the following condition (Dyson & Williams 1980; 
WC89), 

2in^kTH+ = UH^imH^cr'^ + kTk) = nH^mH^Sv^^t/ (8/ti2), (3) 



Note that a uniform density is assumed. 
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where is the molecular gas kinetic temperature and Svtot — {Sln2y/'^{al + -i^^T^Y/'^ 
is the FWHM of an optically-thin molecular line tracing the high density gas around the 
UC HIIs. Svtot includes both non-thermal and thermal components. If S**, the flux of 
ionizing photons, remains largely unchanged during the expansion and the HII region is 
ionization-bounded, one has 



5,(1 - x) = = ^-KRlnlf^^. (4) 



Equations (3) and (4) give 



2 

tot 



So, turbulent pressure can significantly reduce the final equilibrium size of the UC HII 
region. If = 2 km s'^ and ~ 100 K, we have R = 9.5A^^/^Rs (versus R = 5A.3Rs 
with thermal pressure alone). For = 8 x 10^ cm~^ and 5** = 10^^ (for 06 
star), we have Rs = 0.005(1 - xY^^ pc and EMq = 2.5 x 10^°(1 - xY^^ cm-^pc. Thus 
EM = {^YEMq = 2.8 X 108^-^/^(1 - xY^^ cm-^pc and R = 0.05^^/3(^1 - xY^^ pc For 
reasonable dust absorption coefficient x and factor ^, both the predicted size and emission 
measure are reasonable in comparison with those inferred from observations (WC89; 
KCW94), especially considering the expected variation in nH2, and Svtot from region 
to region. In particular, it appears that the initial Stromgren sphere expands by a factor 
of ~ 10, say from ~ 0.01 pc to ~ 0.1 pc, before reaching pressure equilibrium with the 
turbulent ambient molecular gas. The time for the expansion can be roughly estimated from 
the sound crossing time, Texpansion ~ ^ ~ where Cs is the sound speed in ionized 

gas taken as 10 km s~^. Since the average age of the UC HIIs observed is close to 10^ yr 
(WC89; KCW94), most UC HIIs might indeed have had enough time to expand to reach 
pressure equilibrium with their ambient molecular gas (cf. Garcia-Segura & Franco 1996). 

Note that this scenario, just as the classic picture for HII regions, predicts that HII 
regions with larger sizes will have a lower average electron density. This is consistent with 
observations by Garay et al (1993b), Churchwell (1993) and KCW94. Such a size-density 
correlation is somewhat less obvious for UC HIIs with cometary and core-halo morphologies 
(KCW94; Churchwell 1993). But it is somewhat difficult to define the average electron 
density for non-spherical UC HIIs and thus the determined densities for these sources in the 
literature are likely subject to relatively large systematic errors (Stan Kurtz 1996, private 
communication). Further careful observational studies on this aspect will be useful. 



3.2. Velocity Dispersion-UC HII Size Relation 
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Equation (5) indicates that the size of UC HIIs ought to be anti-correlated with the 
total gas pressure or velocity dispersion of the ambient molecular gas if the gas pressure is 
indeed responsible for the compactness of UC HIIs. In the case that Rs (i.e., S^^) does not 
vary much from one region to another, Equation (5) predicts 6vtot « D^fJ,^ . However, the 
surveys by WC89, Garay et al (1993b) and KCW94 imply that has a strong dependence 
on R in the sense that larger UCHIIs tend to have a larger ionizing flux 5* (Churchwell & 
Kurtz 1996, private communication) 0. Assuming that the HIIs are ionization-bounded, the 
— R dependence can be derived from the proposed Ue — R relation, Ue = 'n'eo.ipc(oT^)~") 
where Ueo.ipc is the average electron density for UCHIIs with R = O.lpc. Therefore, 
S'*(l — x) = ^'KR^nlf]2 oc R^~'^°', and we have 



° rriH^ O.lpc 

Churchwell (1993) and KCW94 found that a ~ 0.65 provides a reasonable fit to their 
data for a sample of spherical and unresolved UCHIIs, while Garay et al (1993b) preferred 
a ~ 0.98 for a larger sample of HIIs with various morphological types, some of which have 
relatively large sizes. 

While the size of the UC HIIs can be determined easily from the VLA maps, it is 
not entirely trivial to obtain 5vtot observationally. Specifically, the following criteria must 
be considered in choosing molecular lines for a reliable determination of 5vtot- First, 
the molecular line used must have a relatively high critical excitation density in order 
to sample the dense gas in the immediate neighborhood of UC HIIs. Second, to reduce 
foreground and background confusion as well as line broadening by optical depth effects, 
the molecular line must be optically thin. Third, the spatial and spectral resolutions must 
be adequate. Existing systematic surveys of molecular gas associated with UC HIIs using 
lower transitions of ammonia and CO and isotopes can be immediately ruled out by the 
above criteria. At this time, it appears that the C^'^S survey of 8 UC HIIs by CWKC 
with the IRAM 30m telescope is the only systematic data set in the literature from which 
5vtot can be determined reasonably well for a non-biased sample of UC HIIs. Among the 3 
transitions of C^^S observed, the J = 2 — 1 data with a spatial resolution of 25" is preferred 
based on considerations of its optical thinness, superior velocity resolution and reliable 
detection in all cases. Higher transitions show significantly larger linewidths, indicating 
possible additional broadening due to other physical processes such as optical depth effects, 
and outflows. 



^If S"* does not vary much from region to region, one expects to see a steep rig — i? relation, 
In reality, such steep n^ — R dependence is not seen (Churchwell 1993; Garay et al 1993b; KCW94) 
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Figure 1 shows the data points for all 8 UC HIIs observed in C^'^S J = 2 — 1 by 
CWKC in comparison with the expected correlation for pressure-confined UC HII regions. 
The dotted, solid and dashed lines represent Equation (6) with a = 0.65 (Churchwell 
1993; KCW94) for ^^^^^ = 2.55 x 10-2, ^^^^^ = 4.04 x 10"^ and ^^^^^ = 8.08 x lO'^, 

respectively, i.e. Svtot = (3.16,2.24, 1.78)0^0^^^, where Due = 2R. Despit c a large scatter, 
an anti-correlation between 6vtot and Due is evident and in excellent agreement with the 
theoretical prediction with a = 0.65. In fact, a least-squares-fit to the 8 data points gives 
Svtot{km s-^) = (2.45 ± 1.17) D^^'^^^^-'^^ with a correlation coefficient r = 0.87, which 
is almost identical to the sohd line. Equation (6) with a — 0.98 (Garay et al 1993b) 
does not fit the 8 data points as well, but it is worth noting again that the Garay et al 
sample of UC HIIs contains not only larger HIIs but also non-spherical HIIs for which 
the determined mean electron density is subject to larger uncertainties. Since Svtot does 
not show any correlation with the distances to these Hll regions, the possibility that the 
above anti-correlation results from the different physical sizes probed by the same telescope 
beam at different distances appears unlikely. From Churchwell (1993) and KCW94, we 
have UeQ.ipc ~ 10^ cm^^. Therefore the solid line in Figure 1 implies an average molecular 
gas density n^/j ~ 2.5^ x 10^ cm~^. Given the uncertainties involved, this implied mean 
molecular density is in reasonable agreement with the densities derived for these regions 
(CWKC) and other massive star forming cores (Plume et al 1996), especially considering 
the fact that some UC HIIs may not reside in the densest molecular gas in massive star 
forming cores (cf. Churchwell 1993; Hofner et al 1996; Xie et al 1996). Considering the large 
uncertainty in Due due to the highly uncertain distances to these regions, the consistency 
between the observations and the theoretical prediction for 5vtot — R relation with o; = 0.65 
is surprisingly good. 

Finally, we note that Svtot may include systematic motions such as collapse and outflow 
or expansion, as observations indeed suggest (cf. Ho & Haschick 1986; Peng 1995; Hofner et 
al 1996; Shepherd & Churchwell 1996; Xie et al 1996). The infall of gas provides additional 
confining pressure, just like the turbulent motion. In the extreme case that 5vtot is largely 
dominated by outflowing motion, however, there appears an alternative interpretation for 
the 5vtot — Due anti-correlation in the sense that younger and thus smaller UC HIIs may be 
associated with larger outflowing motions. We feel that this interpretation is less attractive 
given the interferometric observations which indicate that an UC HII is often just one of a 
few massive stars in a cluster on a considerably larger scale, and the energy input due to an 
UC HII alone is unlikely to be dominant (cf. Wilner, Welch & Forster 1995; Xie et al 1996). 



4. DISCUSSION 
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We also note that it would be naive to expect turbulent gas pressure to be the 
dominant confining mechanism for every UC HII region. A desirable aspect of the simple 
turbulent pressure confinement idea is that it does not deny the possible role that any other 
mechanisms may play in addition to the turbulent pressure. 

One natural prediction of the turbulent pressure confinement of UC HIIs is that UC 
HIIs in relatively quiescent molecular clouds must be rare because they will be short-lived. 
Unfortunately, it is difficult to check this prediction observationally, because it has long been 
known that massive stars tend to form in massive molecular cloud cores where turbulence 
is known to be stronger. However, since even an initially quiescent molecular cloud would 
become turbulent as soon as the first generation of stars form and inject kinetic energy 
as outfiows, one interesting speculation is that the first massive stars may develop larger 
HII regions around them and trigger the formation of a later generation of stars due to the 
enhanced gas pressure. The stars of later generations formed in the pressure environment 
would have to spend considerably longer time embedded in their ambient molecular gas and 
dust. 

Systematic surveys of the molecular gas around a large number of UC HIIs using 
optically-thin, high density-sensitive molecular lines with resolution on the order of ~ 10" 
would be highly useful in confirming or rejecting the reality of the proposed velocity 
dispersion-UCHII size relationship and in clarifying the role that gas pressure plays in the 
development of HII regions around newly-formed massive stars. More detailed studies of 
ne — R oi — R relation for UCHIIs of various morphologies and sizes will also be very 
helpful. In particular, it would be highly desirable if high quality multi-transitional data 
could be taken for a large number of UC HIIs with a wide range of sizes, which could then 
be used to determine the density of the molecular gas as well and thus to check if the UC 
HIIs are indeed in rough pressure equilibrium with ambient molecular gas. 
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for their kind encouragement. He further thanks Neal Evans, Pepe Franco, Paul Goldsmith, 
David HoUenbach, Stan Kurtz, Yuan Peng, Debra Shepherd, Prank Shu, Prank Wilkin 
and Qizhou Zhang for useful discussions and Pepe Pranco, Paul Goldsmith, Paul Ho, Luis 
Rodriguez and the anonymous referee for reading the manuscript with helpful comments. 
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Figure Caption 

Figure 1 shows the theoretically expected linewidth-size relationship versus the 
C^'^S J — 2 — 1 data from the survey by Cesaroni et al (1991). Despite a large scatter, the 
data are remarkably consistent with the theoretical predictions for reasonable parameters, 
as discussed in the text. 



